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We study the susceptible- infected model with power-law waiting time distributions P(r) ~ t~", 
as a model of spreading dynamics under heterogeneous human activity patterns. We found that 
the average number of new infections n{t) at time t decays as a power law in the long time limit, 
n(t) ~ f"'', leading to extremely slow prevalence decay. We also found that the exponent in the 
spreading dynamics, /3, is related to that in the waiting time distribution, a, in a way depending on 
the interactions between agents but is insensitive to the network topology. These observations are 
well supported by both the theoretical predictions and the long prevalence decay time in real social 
spreading phenomena. Our results unify individual activity patterns with macroscopic collective 
dynamics at the network level. 



I. INTRODUCTION 



Spreading dynamics represents many real social phe- 
nomena, such as emerging epidemics and information 
flows, motivating its study for many years. Traditional 
epidemic modeling studies have assumed homogeneity in 
the timing of events and in the connectivity patterns of 
contacts [l[. Recently modifications to classical models 
of spreading dynamics have been considered, to account 
for the complex connectivity patterns of social networks 
and the non-Poissonian nature of human activity pat- 
terns 0-01 ■ There have been intensive studies on the 
spreading dynamics in the presence of connectivity het- 
erogeneity, in particular in networks with scale-free de- 
gree distribution P{k) ^ k^'^ with 2 < 7 < 3, that better 
represent the world wide web, co-authorship networks, 
and many other social networks On the other 

hand, the impact of time heterogeneity in human dy- 
namics started to be recognized just recently. Various hu- 
man activity patterns, ranging from tele-communications 
such as E-mail 0, [3|, surface mail web browsing 



12|, instant messaging [13|, and mobile phone 
14 1 to physical contacts probed by wireless de- 
I5L ITgI I or sexual contact survey [13 , have revealed 



calls 
vices 

high heterogeneity regarding the timing of events, which 
is often described by a power-law waiting time distribu- 
tion P{t) ^ T~", where r is the time interval between 
two consecutive activities, a strong departure from tra- 
ditional modeling assumption A few recent studies 
incorporating empirical heavy-tailed waiting time distri- 
bution in the spreading dynamics have demonstrated a 
significantly longer prevalence decay time than expected 
from the Poisson process [l^ [3] ■ 

On the individual level, priority-queue models mim- 
icking human decision process were introduced to repro- 
duce power-law waiting time distributions 0, 0, • In 
these models, prioritization of tasks is considered a major 
cause of heavy-tailed nature of human dynamics. These 
models were extended to include interactions between 
two or more agents inevitable in social system [20l - l22| . 



While these models are successful in reproducing power- 
law waiting time distributions, it is still veiled how the 
activity modeling at the individual level translates into 
the collective phenomena at the population or network 
level. So, here we study spreading dynamics as a rep- 
resentative example of macroscopic phenomena in social 
system using the susceptible-infected (SI) model under 
heterogeneous activity patterns. 



II. SI MODEL WITH POWER-LAW WAITING 
TIME DISTRIBUTIONS 



A. General theory 

To start understanding the impact of power-law wait- 
ing time distribution, we first turn to a general theory of 
irreversible spreading processes in a social network. We 
assume the network has a tree-like structure and charac- 
terize the timing of infection transmissions by the genera- 
tion time distribution g{A), where the generation time A 
is defined as the time interval between the infection of an 
agent in the social network (primary case) and the trans- 
mission of the infection to a neighbor agent (secondary 
case). In this case, an outbreak starting from a single in- 
fected individual ait = results in the average number of 
new infections n{t) = X^d^i Zdg*'^{t) at time t [2^, where 
Zd is the average number of individuals d contacts away 
from the first infected node, D is the maximum of d, and 
g*'^{t) is the d-th order convolution of (/(A) [g*^{t) = g{t) 
and g*'^{t) = dt' g{t')g*'^-'^{t - t')], representing the 
probability density function of the sum of d generation 
times. For the cases where 5(A) ~ A"*^ with 1 < /? < 2, 
in the limit d ^ 1 one obtains g*'^{t) ^ Cp-i{t/td) /td, 

where td ~ Aod^^, Aq is some characteristic time scale 
and Cf^ (x) represents the stable (Levy) distribution with 
exponent ^ concentrated in [0,oo), characterized by the 
asymptotic behavior C^{x) ^ 2;^^^+^) for a; » 1 [2^ . 
From this it follows that, independently of the network 
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FIG. 1. (Color online) (a) The average number of new infections n(t) of the SI model with uncorrelated identical power-law 
P{t) with exponent a = 2.1 (□), 2.5 (o), 2.9 (A), and with an exponential P(r) (v), respectively. All P{t) have the same mean 
waiting time (r) = 2. It decays with the exponent ~ q — 1 for power-law P{r), whereas it decays exponentially for exponential 
P{j). (b) The generation time distribution g(A) for power-law P{t). It decays as a power law with the same exponent as n{t). 
Simulations were performed on a scale-free network with 7 = 3 and size iV = 10*, and averaged over 10* independent runs. 
Initial infected node is selected to be the hub. In both panels, dotted lines have slopes —1.1, —1.5, and —1.9, drawn as a guide. 
Data are binned logarithmically. 



structure, ^ t ^ when t — > oo and 

n{t)^t-f'. (1) 

Therefore, if the generation time distribution decays as a 
povifer law with exponent /3, then the long time dynamics 
of the spreading process will be characterized by the same 
power-law decay. In general, when the generation time 
distribution is heavy-tailed but not power-law, we expect 
that the long-time spreading dynamics decay with the 
same tail as the generation time distribution by the above 
theory. 

B. SI model with uncorrelated power-law activity 
patterns 

The generation time distribution, in turn, can be de- 
termined from the waiting time distribution and the type 
of interactions between the social agents. As the simplest 
scenario of spreading process, let us first consider the case 
when the acitivities of individuals in the network are com- 
pletely uncorrelated and so are the timings of consecutive 
interactions. In such a case, the time series of interac- 
tions can be modeled by a simple renewal process [l^l 
with waiting time distribution P{t). Specifically, let us 
consider the scenario where the agent A is infected at 
time tA, B, who is connected to A, still susceptible at tA, 
and B is infected at time ts after interaction with agent 
A. From the perspective of B, the time A was infected is 
random and therefore, the generation time ts — tA is the 
residual waiting time, the time interval between a ran- 
domly selected time and the time A and B will perform 
their next interaction. 

For the waiting time distribution with finite mean 
(r), the residual waiting time probability density func- 
tion is related to the waiting time probability density 



function, leading to g„„corr(A) = P{T)dT [H. 

Therefore, for the uncorrelated activity patterns with 
P(r) - r"", 2 < a < 3, we have .g„„corr(A) ~ A-^^-^' 
and Uuncorrit) ^ t~^°'~^\ On the other hand, when 

1 < a < 2, P{t) does not have a finite mean so the 
generation time distribution cannot be obtained in the 
same way. In this case the residual time distribution is 
non stationary and, in the long time limit, it approaches 
the limit distribution guncorr{^,t) — > i^/t)Pa{^/t) j23 |. 
where < A < t and Pa{x) = sin7ra/[7ra;"(l + x)]. Fur- 
thermore, in this case the convolutions of (7„„corr(A, 
areoftheformg*5^^„^^(A,t) {I / t)p*^ {l\ / 1) . Substitut- 
ing the latter result into n{t) ~ X^dLi we obtain, 
in the long time limit, n{t) w X^d^i ^dP^i^)- Putting 
all together we thus obtain guncorri^) ~ A"^""'^'"''- for 

2 < a < 3 and nuncorr{t) ^ t-f^^^iorr f-Qj. l < < 3 with 

Puncorr ~ \^a-l, 2<a<3 ^ ' 

To test the predicted power- law decay of n{t), we per- 
form the following numerical simulations. Initially all 
agents are susceptible except for a single infected node. 
Then infected agents infect connected susceptible agents 
following a power-law waiting time distribution with ex- 
ponent a until all agents are infected. For 2 < a < 3, 
n{t) is found to decay as a power law with the expo- 
nent a — 1 for large t as predicted by the theory (Fig. 1). 
For comparison, we also performed the same SI dynamics 
with exponential (Poisson-type) waiting time distribu- 
tion with the same mean waiting time (r) = 2 [Fig. 1(a), 
V]. In this case, n{t) decays exponentially fast, in stark 
contrast to the power-law cases. Specifically, the effec- 
tive duration T of the epidemic process, given by the 
expected infection time of an individual after the out- 
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FIG. 2. (Color online) Schematic illustration of the SI model 
on the PQN. Filled (open) nodes denote infected (suscepti- 
ble). The arrows from an infected to a susceptible node indi- 
cate potential subsequent infection channels. 



break, T = X^t^o *'^(^)/-^' measured to be T^^p ~ 5 
for the exponential P{t), whereas it can become orders 
of magnitude longer as T w 4 x 10^ for the power-law 
P(t) with a = 2.1. Therefore, the power-law waiting 
time distribution indeed impacts the long-time dynamics 
of the spreading process significantly. 



III. SI MODEL ON THE PQN 

In a more realistic scenario, the multi-agent dynamics 
in a social network based on human decision process can 
be modeled by the priority-queue network (PQN) with 
interactions (Fig. 2). We consider the SI model on PQN 
with OR-typc interactions, which can model the spread- 
ing of infectious entity through mutual communications 
or contacts between individuals that can be initiated pri- 
marily by any one of the individuals, such as talking over 
the phone or face-to-face encounter. Specific procedures 
of numerical simulations are as follows. Each node in a 
network is a queue with two kinds of tasks, interacting 
(/) and non- interacting (O) task [lO, HH. A node i with 
degree ki is a queue with fixed length ki + \, having ki 
I tasks and one O task. The / task of the node i paired 
with the node j through the link between them is de- 
noted by lij and the O task of the node i is denoted by 
Oi (Fig. 2). We start with ah nodes being susceptible, 
except for one infected node selected to be the hub. Ini- 
tially a random priority value from uniform distribution 
is assigned to each task. Each step, we choose a node i 
randomly and the highest priority task in the queue is 



executed. At this step, if the highest priority task in the 
queue is Oi task, then it is executed alone. However if 
it is lij task, then not only lij task but also its conju- 
gate task Iji are simultaneously executed. If an infected 
and a susceptible node perform the / task between them, 
the "disease" is transmitted and the susceptible node be- 
comes infected. Then all executed tasks arc replaced by 
new tasks each with a random priority. This infection 
process proceeds until all nodes are infected. We have 
checked that different choices of the initial infected node 
other than the hub do not affect the long-time dynamics 
of spreading processes. 

We first consider the PQN with the scale-free degree 
distribution with 7 = 3. In the long time limit, n{t) is 
found to decay as a power law with exponent approxi- 
mately 2 [Fig. 3(a), o]. The power-law decay of n{t) is 
much slower than the exponential decay in the SI model 
on the random execution queue network [Fig. 3(a), V]. 
Indeed such a longer prevalence time than predicted by 
the Poisson process was observed in real spreading dy- 
namics controlled by human activity patterns, such as 
an E-mail virus outbreak d, [l^ . 

To understand the n{t) ~ t~'^ behavior observed for 
the PQN, we turn to the waiting time distribution of 
PQN, through which the generation time distribution can 
be obtained. It has been shown that there is a local vari- 
ation of the waiting time distribution exponent a in the 
PQN [2l|, which renders the identification of network- 
level exponent apQN nontrivial. Moreover, the highly 
coupled and non-Markovian nature of PQN dynamics 
makes its analytic understanding currently a formidable 
challenge, even for the simplest case with N = 2 (20l.[2l|. 
Lacking exact solutions, however, we could infer its key 
long-time dynamic characteristics from the following sta- 
tistical reasoning: We first observe that different / tasks 
(links) manifest different exponents Pfcfc' {t) ^ t""''''' , 
parametrized by the degrees k and k' of the associated 
nodes. For example, for the / tasks between the hub 
and a node with fc = 1, otkf^^^^.i = 1, whereas a2,i ~ 3 
[Fig. 3(b)]. Putting the contribution of all tasks together 
we obtain Ppqn{t) = J2kk' (^kk/Pkk' (t), where Ckki is the 
fraction of links connecting nodes with degree k and k' . 
The smallest the exponent akk' the biggest is its contri- 
bution to the tail of Ppqm{t), except for the fact that 
links with akk' < 2 cannot contribute statistically be- 
cause of its strong peak of Pkk'{T) at r 1. The net 
result is that links with akk' ~ 2 contribute the most, 
resulting in [Fig. 3(b), o]. 



PQN 



-apQN 



CUPQN = 2. 



(3) 



To proceed with relating the generation time distribu- 
tion with the waiting time distribution, we note that the 
interaction timings between connected pairs are not com- 
pletely independent in the PQN. To see this explicitly, let 
us consider three linearly connected nodes (B — A — C) in 
the PQN. While node A interacts with node B {Iab task 
is being executed), node A cannot interact with node C 
{I AC task has to wait). So the waiting time of Iac be- 
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FIG. 3. (Color online) (a) The average number of new infection n{t) of the SI model on the PQN (o) and on the random 
execution queue network (v). After the initial increase, n{t) after the peak decays with a power-law tail with exponent /3 « 2 
for the PQN (o) and exponentially for the random execution case (V). (b) The waiting time distribution of the whole network, 
Ppqn{t), and the waiting time distributions of links connecting nodes with degree k and k' , Pkk'ij), of the same process for 
the PQN. Both Ppqn{t) and Pkk'iT) decay as a power law, yet with different exponents dependent on the local topological 
position, apQN ~ 2, aige.i ~ 1, 02,1 ~ 3, and 04,4 ~ 2. (c) The generation time distribution (;pqjv(A) of the same process for 
the PQN. It decays as a power law with the same exponent as npQjv(t). In all panels, dotted lines have slope —2, drawn as a 
guide. Simulations were performed on a scale-free network with 7 = 3 and size N = 10*, and averaged over 10^ different initial 
conditions. Data are binned logarithmically. 
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consistent with our numerical findings. 

Finally, we examine the influence of the network struc- 
ture on the spreading dynamics by the SI model over the 
PQN on Erdos-Renyi network and Cayley tree. On both 
network structures, n{t) is found to decay as a power law 
also with /3 = a, n{t) ^ (Fig. 4). Since every node 
except the root in Cayley tree is topologically indistin- 
guishable, all / tasks show the same P(t), which decays 
as a power law with a « 2 for degree fc > 3. These results 
suggest that the network structure exerts minor influence 
on the spreading dynamics in the long time limit and the 
heterogeneous activity pattern is the major determinant 
of the long time dynamics. 



FIG. 4. (Color online) n{t) (left scale) and g(A) (right scale) 
of the SI model on the Erdos-Renyi (ER) network with mean 
degree (fc) = 4 and size = 10* and Cayley tree with branch- 
ing number fc = 4 with 6 generations (N = 1,456). In each 
case, n{t) and g(A) show power-law tails with the same ex- 
ponent. Dotted line has slope —2, drawn as a guide. Data 
are binned logarithmically. 



gins at the execution of Iab task. If the node B has 
been infected, the generation time of infection of C from 
A is the same as the waiting time. Therefore, in the 
PQN the generation time distribution gpQN{A) has the 
same power-law exponent as the waiting time distribu- 
tion Ppqn{t) [Fig. 3(c)], which in turn, according to our 
general theory, leads to npQjs[(t) decaying with the same 
power law as Ppqis[{t), that is 

Ppqn = apQN = 2, (4) 



IV. SUMMARY AND DISCUSSION 

To conclude, we studied the SI model following hetero- 
geneous activity patterns in time. The analytic predic- 
tion indicates that the power-law waiting time distribu- 
tion P{t) leads to the power-law decaying new infection 
number n(t) in the long time limit, with the exponent de- 
termined through the generation time distribution 5(A). 
These theoretical predictions are in good agreement with 
numerical simulation results for two different scenarios of 
spreading process by the SI model on the priority-queue 
networks and that under uncorrelated power-law wait- 
ing time distributions, respectively. In this study, we 
have demonstrated how the macroscopic phenomenon of 
extremely long prevalence time in spreading dynamics 
emerges from individual-level activity model with inter- 
actions. Our results indicate that the heterogeneous in- 
dividual activity patterns significantly impact social dy- 
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namics and therefore it is an essential faetor in modeling 
both the individual and the network level dynamics, such 
as the spreading of human infectious diseases [H and that 
of mobile phone viruses [2^, which is expected to be a 
major problem in the near future. 
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